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A combination of molecular-dynamics (MD) computer simulation and mode-coupling theory 
(MCT) is used to elucidate the structure-dynamics relation in sodium-silicate melts (NSa;) of vary- 
ing sodium concentration. Using only the partial static structure factors from the MD as an input, 
MCT reproduces the large separation in relaxation time scales of the sodium and the silicon/oxygen 
components. This confirms the idea of sodium diffusion channels which are reflected by a prepeak 
in the static structure factors around 0.95 A -1 , and shows that it is possible to explain the fast 
sodium-ion dynamics peculiar to these mixtures using a microscopic theory. 
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Glass-forming mixtures of Si0 2 with an alkali oxide 
such as Li 2 0, Na 2 0, or K 2 are of general interest for 
the study of transport mechanisms in amorphous con- 
densed matter. They all show alkali ion mobility that 
is much higher, often by orders of magnitude, than that 
of the silicon and oxygen atoms forming a tetrahcdral 
network structure. More than 20 years ago, the idea of 
preferential diffusion pathways in the Si-0 network was 
proposed as an explanation for the fast transport of alkali 
ions 0,0 ■ A recent study of sodium silicate melts using a 
combination of molecular dynamics (MD) computer sim- 
ulations and inelastic neutron scattering Q has shown 
that indeed percolating sodium-rich channels are formed 
in the intermediate-range static structure. The sodium 
ions perform a hopping motion along those channels 0, 
and simulation studies on different alkali silicates have 
identified various properties of the characteristic alkali 
ion sites involved in this motion H, 0, HI ■ 

If the presence of intermediate-range order in the static 
structure is intimately related to the high mobility of 
the alkali ions, it should be possible to predict the dy- 
namics, in particular the time-scale separation between 
the sodium transport and the silicon/oxygen "matrix" 
dynamics, from this static structure. Such a structure- 
dynamics relation is, however, not easily formulated in 
terms of a microscopic theory. Here, we discuss a theo- 
retical framework in which this is possible. 

The mode-coupling theory of the glass transition 
(MCT) 0, has been formulated to provide a first- 
principles description of slow dynamics in glass-forming 
liquids. It needs as input the static equilibrium in- 
formation about a system, in its most common form 
only the partial static structure factors S a p{q). The 
MCT equations then yield dynamic correlation functions, 
in particular the normalized dynamic structure factors, 
Fap(q,t) = S a p(q,t)/ 'S a p(q) (also called the coherent 
intermediate scattering functions), and the incoherent 
self-intermediate scattering functions of a given species 
a, F^(q,t). They can be directly compared to scatter- 



ing experiments or simulation results and determine, via 
Green-Kubo relations, the transport coefficients. 

The programme of predicting, using MCT, dynami- 
cal features from the equilibrium structure of the liq- 
uid and testing the predictions against MD simulations, 
has proven successful for some paradigmatic dense liquids 
like the hard-sphere system |l0j , binary mixtures of hard 
sphere s ITDl , Lennard- Jones particles [13, UM > a model of 
NiZr4 Ibi and models of more complicated dense glass 
formers~ fl5l flil Il7j . In these systems, the agreement was 
qualitative and often even quantitative. For a computer- 
modeled silica melt, MCT was able to predict the form 
factor of the glass, i.e. the height of the plateau in the dy- 
namic structure factors, qualitatively, and even to within 
a few percent if in addition to S n p(q) static triplet cor- 
relations were taken into account |l8| . 

In this work, the MCT equations are solved for sodium 
silicates, using MD-simulated static structure factors as 
an input. The resulting dynamic structure factors are 
compared to the ones obtained from MD. We demon- 
strate that, at least on a qualitative level, the theoretical 
framework of MCT allows to describe the large time- 
scale separation in sodium silicates. This establishes a 
structure-dynamics relation based on static pair correla- 
tions. It confirms that the high mobility of sodium ions 
in sodium silicates is indeed intimately related the pres- 
ence of intermediate range order (i.e. the aforementioned 
channel structure) as reflected by a prepeak in partial 
static structure factors. 

MD simulations were performed for sodium silicates 
with different compositions, (Na 2 0)(Si0 2 ) x [NSa;], with 
x = 2, 3, 5, 20, and 40. The simulations were all done 
at constant density p = 2.37 g/cm 3 , using N = 8064 par- 
ticles and a slight modification of the pair potential by 
Kramer et al. |l9|, which is based on ab initio calcula- 
tions. For details, we refer to Ref. [13, Based on 
the thus obtained S a p(q) for NS2 and NS20, the MCT 
equations of motion for S n fl (q, t) were solved numeri- 
cally as outlined in Ref. [23j , with discrete wave numbers 
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FIG. 1: Partial structure factors for the indicated sodium 
silicates at T = 3000 K, obtained from MD simulations, a) 

S'NaNa(c7), t>) SsiNa(?)- 

0.1 A _1 < q < 20 A _1 of step size Aq = 0.1 A~ , ap- 
plying a five-point running average to the input. The 
MCT solutions are fully determined only through S a /3(q) 

and the static triplet correlation function c a p 7 (q, k). We 
set the latter to zero thus approximating the three-point 
correlations by a convolution approximation. 

The structure factor input is exemplified in Fig. ^ 
where we plot SNaNa(<?) and S'siNa(f?) for different sodium 
concentrations at the temperature T = 3000 K. Note 
that even at this high temperature, typical "matrix" re- 
laxation times are already in the range of 10-100 ps. The 
channel structure formed by the alkali ions on interme- 
diate length scales, at about 6-7 A, is clearly visible as 
a prcpeak around q = 0.95 A . This prepeak becomes 
more intense with increasing sodium concentration, and, 
remarkably, its location does not change significantly over 
the whole considered sodium concentration range. This 
is in agreement with our recent finding that the prepeak 
reflects a characteristic intermediate-range chemical or- 
dering of sodium that leads to the formation of sodium- 
diffusion channels. However, in Fig. ^ no obvious scale 
separation hints towards the peculiar sodium dynamics. 
In MCT for single-component systems, relaxation times 
are typically correlated with peaks in S(q); thus, from a 
peak in the diagonal clement S aa (q) , one might expect a 
slowing down of the dynamics around the prcpeak posi- 
tion, but nothing particular for the self-relaxation times. 

In the comparison of dynamical quantities between MD 
and MCT, one has to note that the theory predicts a di- 
vergence of relaxation times at T = T c : the idealized 
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FIG. 2: Normalized sodium-sodium intermediate scattering 
functions, JNaNa(g) (left panel), and sodium-self intermediate 
scattering functions, F-§) a (q,t) (right panel), for sodium disili- 
cate (NS2). The symbols are MD results for T = 2100 K (from 
Ref. II); the solid lines are MCT results for Tmct = 3410 K; 
q values are indicated. 



glass transition, which is smeared out in experiment. A 
scaling analysis of the MD data gives a value for T c , 
which is lower than the value calculated within the the- 
ory. Hence, in order to compare dynamical data that 
show the same degree of slowing down, i.e. the same sep- 
aration between the structural relaxation time scale and 
the time scale of molecular short-time motion, we com- 
pare MD data for temperature T with MCT data at a 
different (higher) temperature Tmct- This is then the 
only free parameter in our comparison. 

Note that, despite the quantitative error in the value 
of T c , the qualitative change of T c with varying sodium 
concentration is captured by MCT: from a scaling analy- 
sis of the simulation data, one gets T c ps 3330 K for silica 
Hi] and 2000 K for NS2 [13]. The respective values we 
calculate within the theory are T c ps 3983 K and 3105 K. 

The large difference in time scales between the ma- 
trix relaxation and that of sodium transport, on almost 
all length scales, is best evidenced in the sodium-sodium 
(coherent) and the sodium-self (incoherent) correlation 
functions H, F aa (q,t) = S aa (q,t) / S aa (q) and F°(q,t) 
in the case a = Na. The former follows the slow matrix 
relaxation, while the latter decays much faster, giving rise 
to fast sodium diffusion Q. MD results at T = 2100 K 
together with MCT results at T M ct = 3410 K for the two 
quantities at selected wave numbers q are shown in Fig.0 
We have chosen the latter temperatures for the compar- 
ison, since at these temperatures the coherent relaxation 
times as obtained from MD and MCT are very similar 
at q = 0.95 A (see also below). At this wave num- 
ber, the time-scale separation is most pronounced: while 
the sodium-sodium correlations decay at t ps 10 4 ps, the 
sodium-self correlations decay at t s ps 10 2 ps, i.e. two or- 
ders of magnitude faster. MCT underestimates this sep- 
aration, but it still gives a factor of about 50 between the 
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FIG. 3: Structural relaxation times r a (q) for coherent den- 
sity relaxation functions F aa (q,t) and rj a (g) for the incoher- 
ent density relaxation function F^ a (q,t), a) To(q) and tq((?), 
b) iNa(g) and TN a (g). As indicated, symbols are simulation 
results (T = 2100 K) and lines are MCT results (T = 3410 K). 
See text for the definition of r a (q) and r^iq)- 



two time scales. Note that in the hard-sphere system, the 
coherent and incoherent relaxation times agree within a 
factor of 4 at comparable intermediate q. Hence, the the- 
ory captures the peculiarity of the fast sodium dynamics 
in NS2 qualitatively correctly. At the higher wave num- 
ber q — 2.05 A , the time-scale separation is smaller, 
but still t/t s w 20 for the MD data. The MCT result, 
t/t s w 14, here is even closer to the simulation value. 

Overall, the theory reproduces well the qualitative 
facts seen in the MD simulations. There are two ma- 
jor types of discrepancy visible: first, the heights of the 
plateaus seen in the correlation functions at intermediate 
times differ. They are connected to the glass form fac- 
tors, for which in the case of silica [lij the agreement 
was quantitative with static triplet correlations taken 
into account, but qualitative only without them (as in 
our case). We therefore expect this first problem to be 
curable by providing better static structure input. The 
second problem is the shape of the correlators, especially 
the sodium-self correlators at long times: for example, for 
q = 0.95 A , MCT predicts a long, extremely stretched 
decay from 0.6 to zero over more than 2 decades in time, 
which is not observed in the simulation. 

To test the wave-vector dependence of the dynamics 
in more detail, we turn to a discussion of structural 
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FIG. 4: Normalized sodium-sodium intermediate scattering 
functions, JNaNa(g) (left panel), and sodium-self intermediate 
scattering functions, F^^(q,t) (right panel), as in Fig. H but 
for NS20. MD results (symbols, from Ref. H|) are for T = 
2500 K, MCT results (solid lines) for Tmct = 4000 K. 



relaxation times for the coherent, T a (q), and the inco- 
herent correlators, r*(q). We define them through fits 
of stretched exponentials to the correlation functions, 
F aa (q,t) « A a (q)exp[-(t/T a (q)) pa ^] at long times, 
with (3 a (q) < 1 and A a (q) < 1, and similarly for F*(q,t). 
This yields a reasonable description of the dynamics for 
t > 1 ps. In Fig. |2| coherent and incoherent relaxation 
times are plotted for oxygen and sodium. MD and MCT 
results arc shown as symbols and lines, respectively. The 
coherent relaxation times are, at fixed q, all of the same 
order of magnitude, and they almost coincide around 
q ?s 1 A , the region corresponding to the main struc- 
tural features of the system (note that this holds also 
for Tgj(g)). The values obtained from MCT and those 
from MD are in good agreement given that we only have 
one adjustable parameter. In particular, their subtle q- 
dependence is reproduced. The oxygen-self relaxation 
time To(<7) is of the same order as the corresponding 
coherent one, To(q), except for q — > 0. This is what 
one expects in a typical glass-forming liquid. On the 
contrary, the sodium-self relaxation time T^ a (g) is much 
smaller than all other structural relaxation times, both 
in the MCT and the MD results. Even if the MCT re- 
sults for this quantity are systematically higher than the 
MD data (reflecting the difference in shape of the sodium 
self-correlators in Fig. |2J) , this confirms that MCT indeed 
predicts the fast-ion dynamics in NS2 from the equilib- 
rium structure alone. 

We now investigate the effect of sodium concentration 
on the dynamics by discussing simulation and MCT re- 
sults for NS20. Fig. 0] shows the coherent and incoher- 
ent sodium-density correlators for this system, similar to 
the NS2 results shown in Fig. □ For NS20, results at 
higher T than in NS2 are shown, in order to compare 
states which show the same degree of slowing down, i.e. 
r Na = O(10 3 ps) for q w 0.95 A" 1 . The main effect of 
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reducing the sodium concentration from NS2 to NS20 
on the dynamics is that the sodium-self relaxation now 
is appreciably slower, as is evident from the right-hand 
panel in Fig. 0] The relaxation time obtained from the 
MD data is T^ a (q = 0.95 A 1 ) ps 10 3 ps, similar to the 
7Na(9) value at that q, i.e. all relaxation times, including 
that of the sodium-self correlations, are now of the same 
order of magnitude at fixed q. MCT again somewhat 
overestimates the sodium-self relaxation times, but the 
qualitative sodium-concentration effect is well explained 
by the theory. 

Also for NS20, the shape of the relaxation curves pre- 
dicted by MCT for the sodium-self correlators disagrees 
with that observed in MD. But here, the discrepancy 
is found largely in the form of the relaxation towards 
the plateau, while the further deviations visible in Fig.0] 
are ascribable to the error in those plateau values as ex- 
plained above. The observation of deviations at relatively 
short times, compared to the structural relaxation time, 
is typical for comparisons between MCT and Newtonian 
dynamics MD simulations fill fl"3| . 

In summary, we have used MCT in conjunction with 
computer-simulated static structure factors for different 
sodium silica melts. The theory is able to explain the de- 
coupling of time scales between the sodium-self dynamics 
and the much slower network relaxation, as well as the 
dependence of this decoupling on the sodium concentra- 
tion. The ratios of relaxation times are in qualitative 
agreement with the MD results. Since MCT only makes 
use of the partial structure factors S a p(q), this shows 
that indeed a structure-dynamics relation holds for the 
fast ion transport in these alkali silicates, and that MCT 
provides a possible microscopic explanation for it. The 
description is not yet on a quantitative level, but one 
can expect that taking into account more detailed infor- 
mation about the equilibrium static structure, i.e. static 
triplet correlations, will improve on at least some of the 
deviations pointed out above. 
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